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Abstract 

It has been known for a long time that the ground state problem of random 

magnets, e.g. random field Ising model (RFIM), can be mapped onto the max- 

flow/min-cut problem of transportation networks. I build on this approach, 

relying on the concept of residual graph, and design an algorithm that I prove 

to be exact for finding all the minimum cuts, i.e. the ground state degeneracy 

of these systems. I demonstrate that this algorithm is also relevant for the 

study of the ground state properties of the dilute Ising antiferromagnet in a 

constant field (DAFF) and interfaces in random bond magnets. 
PACS numbers: 02.70.-c, 75.10Hk, 02.60.Pn, 02.10.Eb 
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I. INTRODUCTION 



The statistical physics of random and frustrated systems has received a considerable 
deal of attention in recent years ffl. The presence of quenched disorder has been found to 
greatly change the bulk and interface properties of a variety of systems as compared with 
their 'pure' counterparts, leading to new and very interesting equilibrium and nonequilibrium 
phenomena. Unfortunately the progress has been many times rather slow, primarily because 
random systems pose sometimes insurmountable theoretical difficulties even to the most 
stubborn theorists. Computer simulations have played and continue to play an important 
role in the field, being at times the only guide through a very complicated energy landscape. 
While the traditional Monte Carlo method |2| proved its usefulness again and again, it 
was soon realized that other approaches should be considered, depending on the nature of 
the problem at hand. Since then a variety of algorithms, previously known only within the 
computer science community, have been successfully brought to bear on numerous statistical 
mechanics problems with quenched disorder, from spin glasses to rigidity percolation 
Q. Such algorithms, generally known as combinatorial optimization algorithms, have been 
typically used to find the exact T = ground states of the system being studied, completely 
avoiding the equilibration problems specific to the Monte Carlo simulations. 

In the following I will focus on a single class of such algorithms, network flow algorithms, 
that have been put in the limelight by the work of Ogielski [||], who applied them to the 
study of the random field Ising model (RFIM). Since then the same method was also suc- 
cessfully applied to the study of equilibrium interfaces in disordered systems |||7| , becoming 
an important tool for the physicists working in the field. The method is generally based 
upon mapping the system being studied onto a network of capacities through which an in- 
compressible fluid obeying local mass conservation flows. The problem of finding the ground 
state turns out to be equivalent to finding the maximum flow that can be pushed through the 
network between two special nodes, the source and the sink, the so called max-flow/min-cut 



problem of operations research lq-10]. The advantage of this approach is that polynomial 
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time algorithms have been developed for this problem |9HTT||, some of which are much older 
than the field of random systems. 

One of the problems that was usually not addressed using these algorithms and an 
important problem in the physics of random systems is the degeneracy of the T = ground 
states. In terms of the associated network flow problem this is the question of minimum 
cut degeneracy. An approximate algorithm dealing with this issue was proposed in ]12| and 



applied to the RFIM problem. In this paper I build on the max-flow/min-cut approach, 
relying strongly on the concept of residual graph [|ll|], and design an exact algorithm for 



finding all the minimum-cuts (or equivalently all the ground states for a certain class of 
systems) . 

The organization of the paper is the following: for the sake of completeness Section 
II introduces the mapping of the ground state problem to the max-flow/min-cut problem 
along with the network flow terminology and two combinatorial results that will be used 
in the design of the algorithm, Section III describes the degeneracy algorithm, Section IV 
presents a number of applications of the algorithm and the Appendix includes the proofs to 
Propositions 1-4. 

II. GROUND STATES USING MAX-FLOW/MIN-CUT ALGORITHMS 

In the following I will present the mapping of the ground state problem to the max- 
flow/min-cut problem and then proceed to describe what kind of information one can extract 
from this mapping. I will use the T = interface problem in the random bond Ising 
ferromagnet to illustrate the method because it is somewhat easier for the unfamiliar reader 



to understand it intuitively ( see also |13j for a review). 
The Hamiltonian of the system is 

H=-Y,JijVi<7j (1) 

(ij) 

where Oi = ±1 and Jjj > are ferromagnetic couplings between neighboring spins. Jjj 
are fixed independent identically distributed random variables - quenched randomness. If 



the system is a <i-dimensional cube, an interface with dimension d — 1 can be induced by 
using periodic boundary conditions along d — 1 directions and setting the spins in the two 
d — 1 dimensional hyperplanes that represent the boundaries of the system along the last 
direction to +1 and —1, respectively. The interface that will form in the system between 
the +1 and —1 hyperplanes will generally be rough, because it will wander in order to break 
the weakest bonds. The energy is a minimum over all the possible spin configurations (with 
the ±1 boundary spins fixed) and therefore the problem of finding the minimum energy 
configuration(s) would appear to be computationally very hard, even for small system sizes. 
As it turns out, this is not the case if one takes advantage of the similarity between this 
problem and the max-flow/min-cut problem, very well known in the study of transportation 
networks [jSHlQfl. The idea is the following: two new extra sites are introduced, a source node 
s that is connected to all the spins of the +1 hyperplane, and a sink node t connected to all 
the spins of the —1 hyperplane. The ferromagnetic constants coupling the source node a s 
and the sink node a t to their corresponding hyperplane are chosen to be strong enough so 
they are not broken in the ground state. Then by setting a s — +1 and <j t = — 1 an interface 
is induced as before. Now we view the system, including s and t, as a graph whose arcs are 
the bonds between the spins. The arcs have forward and backward capacities equal with the 
corresponding coupling constants, = Cji = Jij, or we can imagine that the nodes i and 
j are connected by both forward and backward arcs with capacities Gy and c^, so this is a 
directed graph. ( The constraint that c i:; - = is dictated by the physics and is not specific 
to the general max-flow/min-cut problem. It can be relaxed for the interface problem, but 
not for the random field problem.) We define the set of nodes as N and a partition (X, Y) 
of the nodes as 

X EE {l G Nfa = +1} (2) 

Y = {i G N\a, = -1} (3) 

Then X U Y = N, I H 7 = 0, sGX and t G Y. The knowledge of such a partition 
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determines the energy of the spin configuration and the position of the interface. This is 
readily seen if we write the Hamiltonian of the system as 

H = - E Jv~ E J ij+ E Jij = H + 2 e J v ( 4 ) 

(iJ)eA(X) (i,j)eA{Y) (i,j)&A{X,Y) (i,j)eA(X,Y) 

where H is the energy of the fully aligned system, H Q = Yu,j) an d we defined A(X) = 
G X,j G X}, A(Y) = G YJ G Y} and A(X,Y) = G X,j G Y}. 

Thus the problem of finding the ground state interface, which has minimum energy, is 
equivalent to finding a partition (X, Y), also called a cut, that minimizes Y(i,j)EA(x,Y) Jij 
- minimum cut. (Note that H is a constant for a given random sample.) If we imagine 
fluid flowing through the network from the source s to the sink t, the minimum cuts are the 
bottlenecks of the network of capacities because they determine the maximum flow that can 
be pushed through the network from s to t. It is useful at this point to put the problem in 
a more general setting. For this a number of definitions are in order. 

A transportation network is a directed graph G C (N, A) with two special nodes, s - source 
and t - sink or target; N is the set of nodes and A the set of arcs. The directed arc (ij) 
going from node i to j has capacity Cy > 0. (For the sake of clarity we assume that if the 
arc (ij) exists and has capacity Cjj > then (ji) also exists and has capacity Cji > 0.) A 
flow through the network G C (N, A) is a set of numbers {fij}, each corresponding to an arc 
in A, subject to the following feasibility constraints: 

< fij < Cij capacity constraint (5) 
ti = E fji ~ E fv = ^ local flow conservation (6) 

{j\(ji)£A} ' {j\(ij)£A} 

for all the nodes j G N — {s, t}, and 

-e s = e t = f (7) 

where / is called the value of the flow. (Note that — e s = et follows from the flow conservation 
0.) The maximum flow problem of network flows is concerned with finding the flow {fij} 
through the transportation network G c that has maximum value /. 
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For any feasible flow through the network we define the residual graph as the graph 
Gf(N,A) with positive residual capacities of the arcs, 

r H = °ij ~ fij + fji > ( 8 ) 

where — fji is the net flow from i to j. (Note that r^ > follows from the capacity 
constraint, and also that it is possible that c^- = and > 0, when Cji > 0.) An 
augmenting path is a directed path from s to t in the residual graph Gf. 

A cut is a partition of the nodes set N into two subsets X and Y, denoted by (X,Y), 
with s G X and t G Y, such that X U Y = N and X H F = 0. The capacity of a cut is 
defined as: 

C (X,r)= £ Cy (9) 

{{ij)eA\iex,jeY} 

In the following I will concentrate on the case of a transportation network G C (N, A) in 
which if the arc (ij) exists and has capacity > than the arc (ji) also exists and has 
capacity Cji > 0. This kind of network is the most relevant one in the physics applications 
that I described before. For these networks the following two propositions can be proven 
(the proofs are contained in the Appendix), that can be used to design an algorithm that 
finds all the minimum cuts in the network G C (N, A). 

Proposition 1: If {fij} is a maximal flow and (X, Y) a minimum cut then ry = for all 
arcs {(ij) G A\i G X, j G Y}. 

Proposition 2: If f m ax > a cut (X, Y) in G c is a minimum cut if and only if it is a 
directed partition in G™^, the residual graph for a maximal flow. 

The calculation of an actual maximal flow through the network can be done using poly- 
nomial time algorithms. The first such algorithm, the augmenting path algorithm |J, was 
proposed originally by Ford and Fulkerson and it is also a way to prove the max-flow/min-cut 
theorem (see Appendix). However, much faster algorithms have been developed in recent 
years, in particular push-relabel methods with global updates, that allow one to deal with 
much bigger systems than before [|K|. These algorithms can be and have been used in such a 
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way as to improve the computational speed by providing only the value of the maximal flow 



and a minimum cut, but not an actual maximal flow through the network JTTJ] . However, 
the knowledge of an actual maximal flow through the network is crucial for finding all the 
minimum cuts. Fortunately, such a calculation can be made without a major loss of speed 

Ga- 



in. A DEGENERACY ALGORITHM 

In the following I will use Propositions 1 and 2 to design an algorithm that finds all 
the minimum cuts. The algorithm will be aimed at finding the set of arcs {(ij) G A\i G 
X,j G Y,for all minimum cuts (X, Y)}, denoted thereafter by A mc , and a procedure for 
determining the actual minimum cuts. 

Let us assume that we constructed a maximal flow {/™ ax } through the network G C (N, A) 
using an appropriate algorithm and let G™ ax be the associated residual graph. We define: 

Y t = {all nodes that can reach the sink along a directed path in G™^} (10) 

X t = N-Y t (11) 

Then (X t , Y t ) is minimum cut that we will call the T cut (this follows from Proposition 2). 
This minimum cut has the property that for any other minimum cut (X, Y), Y H Y t = Y t . 
{Proof. Let us assume Y D Y t = Y\ ^ Y t ; then X\ = Y t — Y\ C X. Because (X,Y) is a 
directed cut (Proposition 2), for any node i G X there is no directed path from i to the sink 
t in G™ a:c , which must also be true for any % G X\. But X\ C Y t , so this is a contradiction.) 
We also define: 

X s = {all nodes that can be reached from the source along a directed path in G , p ax } (12) 

Y S = N-X S (13) 

Then (X S ,Y S ) is also a minimum cut that we call the S cut. This minimum cut has the 
property that for any other minimum cut (X, Y), X fl X s — X s . Now we define 
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Z = N-X s -Y t (14) 

Then if Z = a single minimum cut exists, otherwise there are at least two. We will be 
concerned with the non-trivial case Z ^ 0. 

At this point our knowledge of A mc is summarized in Fig. [I]. The arrows stand for 
possibly more than one arc of G(N, A), and all these arcs are included in A mc . (Note that 
these arcs are saturated by the flow, that is rij = 0, so in G™ ax only the arcs (ji) are 
present). By the construction of Z we also know that all the remaining arcs that make up 
A mc connect nodes that are included in Z. 

The first part of the algorithm is aimed at finding the disconnected pieces (clusters) that 
make up Z, and is therefore called cluster counting. A look at Fig. |], which is an application 
of the algorithm to the bimodal RFIM, clarifies the significance of this step. The idea is 
to determine the connectivity of Z in GJ ,ax using as connectivity rule (r^ ^ or ^ 0) , 
which is equivalent to (cy 7^ or Cji 7^ 0). This procedure does not reveal any new arcs of 
Amc, but shows how the minimum cuts are constructed using the arcs of A mc that we know 
at this moment. The well known Hoshen-Kopelman algorithm which is a an efficient 



cluster labeling procedure, can be readily adapted for the cluster counting task, and allows 
such a calculation to be made in a time proportional with the number of nodes in Z ||i~5 



Fig. |2] summarizes our knowledge of the minimum cuts after this step. It is easy to see that 
the number of minimum cuts that we can construct at this time is 2 n<yZ \ where n(Z) is the 
number of independent clusters making up Z. Our search for the remaining minimum cuts 
is then reduced to finding the minimum cuts in each of the independent clusters. If f ma x = 
the algorithm can be stopped here as everything is known about the minimum cuts. Any 
partition of the n(Z) clusters is a minimum cut and the total number of minimum cuts is 
2"( z ). (Note that if f max > the 'zero' step of the algorithm should be to find the cluster 
percolating from s to t with the above connectivity rule. All the other clusters contribute in 
a trivial manner to the minimum cuts, therefore in the following we assume that they have 
been already discarded.) 
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In the second part of the algorithm we perform a subcluster counting procedure on each 
of the independent clusters found at the first step. The connectivity rule that we use is 
(fij 7^ and Tji 7^ 0). It is clear, using Proposition 1, that the subclusters so obtained 
can only be on one side or the other of a minimum cut, i.e. the set A m c is included in 
the set of arcs that connect these subclusters to one another (we formally call X s and Y t 
subclusters). In order to eliminate the overcounting, we apply a third procedure called 
subcluster coagulation. 

At the end of the second step of the algorithm the subclusters are connected with each 
other possibly through multiple arcs with = and > 0. The subcluster coagulation 
procedure is applied iteratively on each of the independent clusters and consists of the 
coagulation of subclusters that make up a directed cycle. A directed cycle is an ordered 
sequence of subclusters, Si, S 2 , ■ S m , such that for all Sk, k = 1, ...,m, exist ik, jk nodes, 
h, jk £ Sk, not necessarily distinct, with the property that r ik j k+1 = and rj k+lik > 0, 
where S m+ i = S\. The idea is that, according to Proposition 2, subclusters making up 
a such directed cycle cannot be on different parts of a minimum cut, therefore the arcs 
connecting them are not included in A mc (see also the discussion following Proposition 4). 
A Hoshen-Kopelman type algorithm, in which the subcluster coagulation is achieved through 
relabeling of the subclusters, is again rather efficient at handling this task, but managing 
the data structure requires rather intricate coding. 

At the end of the algorithm we have constructed a supergraph like the one in Fig. |3|, 
which we denote by Q{M,A), with single directed arcs and no cycles. The nodes are now 
the subclusters and the arcs stand, as before, for possibly more than one arc of G(N,A). 
Formally, the arc (JJ) going from the subcluster I to subcluster J is defined by: 

(IJ) = {(ij)eA\ieI,je J} (15) 

and by the construction of the algorithm = and > for all arcs (ij) 6 (I J)- Then 
we define A sg = LU(i"J), the set of arcs of G(N, A) represented by the arcs of G(Af, A). The 
following proposition is then true (see Appendix for the proof): 
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Proposition 3: At the end of the algorithm A sg = A m c- 

When the algorithm terminates A mc is known and moreover, the problem of counting 
the minimum cuts is reduced to finding all the directed partitions in a directed, much 
smaller supergraph - G(N , A) - with single arcs and no cycles. Simple enumeration is 
therefore feasible if the independent clusters are not too big. The degeneracy (total number 
of minimum cuts) D can be written as: 



where d(i) is the degeneracy associated with cluster i, d{i) > 2. 

One additional result can be proven (the proof is contained in the Appendix), that was 
also probably known by Ford and Fulkerson ||, that further clarifies the meaning of A mc - 

Proposition 4'- A mc is the set of arcs of G(N, A) that will be saturated by all the maximal 
flows (the 'weak' links of the network). 

The above proposition makes it easier to understand intuitively the significance of the 
subcluster coagulation step of the algorithm. It is clear that the flows between subclusters 
are all saturating and also that the flow conservation holds for each subcluster as a unit. 
Therefore, if for a certain maximal flow we identify a saturating flow cycle between subclus- 
ters, we can reduce the flow around the cycle by its smallest value between two subclusters 
making up the cycle, while keeping the overall flow maximal. This implies that none of 
the arcs making up the cycle is necessarily saturated when the flow is maximal, therefore, 
according to Proposition 4, they are not part of A mc . Furthermore, because of the way the 
subclusters have been constructed, none of the arcs connecting the subclusters making up 
the cycle is part of A mc , so the subclusters can be coagulated. 



Our interest in designing a degeneracy algorithm arose in connection with our desire to 
study the ground state properties of random magnets. We first used the above algorithm 




(16) 



IV. APPLICATIONS 
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to study the ground state structure of the two-dimensional ±h random field Ising model 
(RFIM), that we expected to have a large degeneracy. The Hamiltonian of the RFIM is: 

H RF im — —J^ViVj -^hidi (17) 

(ij) * 

J > and hiS are independent random variables drawn from a symmetric distribution P(hi), 
with (hi) = and (hf)^ = h. In the case of the ±/t RFIM the h\s are, with equal probability, 
+h and —h. This problem is mapped onto the network flow problem by connecting the sites 
with positive fields to the source and the ones with negative fields to the sink through arcs 
with capacity h (see for example |13| for more details). 



The structure of a ±/i RFIM ground state is shown in Fig. |j: two frozen ferromagnetic 
domains of 'up' and 'down' spins and a number of isolated clusters that can be flipped (col- 
ored) independently of each other to generate new ground states. There is also a degeneracy 
associated to flipping certain groups of subclusters inside a cluster. Surprisingly, domains 
that can be flipped without changing the energy exist even if h/J is irrational. In this case 
they have zero field energy and the same exchange energy in the 'up' and 'down' states, and 
are located at the boundary of the frozen domains. The smallest such clusters have only 
two spins ( see Fig. As a result of this structure the ±/i RFIM has a strictly positive 
entropy for a range of h/J, and related with it a new order parameter, the paramagnetic 
response associated with the orientation of the above domains. This result may be relevant 



to universality issues that are currently being debated flq|17 |. 

A similar structure is found for the dilute Ising antiferromagnet in a constant field 
(DAFF), Fig. |5|, which is believed to be the experimental realization of the RFIM: 

Hdaff = J titjOiOj ~ h J2 e i a i ( 18 ) 

(ij) i 

J > and €i are quenched independent random variables, -P(ej) = p5(e« — 1) + (1 — p)5(e i ) ) 
< p < 1. Detailed results are being reported elsewhere [IB] . 



For the sake of clarity one of the clusters in Fig. [| and its subcluster graph representation 
as given by the algorithm are shown in Fig. |^. The connections to X s and Y t are not shown; 
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any directed partition is a minimum cut. Note that the arrows separating such a partition 
point from the 'up' subclusters to the 'down' subclusters. 

We also applied the algorithm to the study of the ground state interfaces in the bond- 
diluted Ising model. The Hamiltonian is the one in Section II, with Jjj being J > with 
probability p and with probability 1 — p. The ground state structure is shown in Fig. [7], 
with neighboring subclusters having different colors. These subclusters can be thought of as 
the excitations of a single interface, as in Fig. [7|, or multiple interfaces as in Fig. || The size 
distribution of these excitations is a power law. Detailed results will be reported elsewhere 

In conclusion, we designed an exact algorithm that finds all the arcs of a flow network 
that are part of a minimum cut and allows the effective construction of all the minimum cuts. 
This algorithm is relevant for the study of the ground state properties of the random field 
Ising model (RFIM), dilute Ising antiferromagnet in a field (DAFF), interface properties 
of certain Ising models with bond or site disorder and for other physical problems that 
can be mapped onto network flow problems and where the ground state is expected to be 
degenerate. 

After this work was completed Bruce Hendrickson brought to my attention the work of 
Ball and Provan [20[|, that contains the idea of constructing an acyclic graph that can be 



viewed as a compact representation of all minimum cuts (see also f2T|). The authors also 
present a polynomial algorithm for counting the directed cuts in a planar acyclic graph. It 
appears that these problems are of continuing interest for the study of network reliability. 
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VI. APPENDIX 



In the following we present the proofs of the Propositions 1-4 that were quoted in the 
main text. They rely in part on the max-flow/min-cut theorem of Ford and Fulkerson ||[|, 
probably the most important result in network flows: 

Theorem (Ford and Fulkerson, 1956): In a transportation network G C (N, A) the max- 
imum value of / over all flows {fij} is equal to the minimum value c(X,Y) over all cuts 
(X,Y). 

Proposition 1: If {fij} is a maximal flow and (X, Y) a minimum cut then r^ = for all 
arcs {(ij) G A\i G X,j G Y}. 

Proof. Let {fij} be a maximal flow with value f m ax and (X,Y) a minimum cut. We 
have / = f max = Y,{(ij)eA\iex,jGY}(fij ~ fji) b Y conservation of the flow [|], and also 
fmax = J2{(ij)eA\iex,j€Y} c ij> by the Ford-Fulkerson theorem. This immediately implies 
Y^{(ij)eA\iex,jeY} r ij = 0> an d because r^ > Proposition 1 follows. 

Proposition 2: If f max > a cut (X, Y) in G c is a minimum cut if and only if it is a 
directed partition in GJ' ax , the residual graph for a maximal flow. 

Definition: A partition (X*, Y*) of the nodes N, X*UY* = X, X* n Y* = 0, is directed 
if the arcs connecting X* and Y* all have the same direction, for example going from Y* to 
X*, i.e. Tij = for all arcs {(ij) G A\i G X*,j G Y*} and 3{j*i*) G A, i* G X*, j* G Y*, 
rj*i* > 0. 

Proof. If (X, Y) is a minimum cut we have from Proposition 1 = for all arcs 
{{ij) G A\i G X, j G y}. Now if we also assume Tji = for all arcs {{ji) G A\i G X, j G F} 
this implies c^- = Cjj = for all arcs {(ij) G A, (ji) G A|i G X, j G F}, and further that 
fmax = from the Ford-Fulkerson theorem, which is a contradiction. Thus the direct impli- 
cation is proved. 

Conversely, let us now assume that (X*,Y*) is a directed partition in GJ - ^, i.e. r^ = 
for all arcs {(ij) G A\i G X*,j G Y*} and 3(j*i*) G A, i* G X*, j* G Y*, r jH * > 0. In 
order to prove that (X*,Y*) is a minimum cut we have to prove first that (X*,Y*) is a 
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cut, i.e. s E X* and t E Y*. Let us assume that this is not true. Then either a)(s E Y*, 
t E X*), or b)(s EX*,tE X*), or c)(s E Y*, t E Y*). We will show that all these lead to a 
contradiction. 

a) (s E Y*, t E X*). In this case (Y*,X*) is a cut and by the conservation of the 
flow J2{(ij) € A\i€X*,jeY*}(fji - fij) = fmax > 0. But r {j = implies fji - f i:j = -c i:j , so 
52{(ij)eA\iex*jeY*}(fj% ~ fij) = ~E{(ij)\iex*,jeY*} < 0, and therefore a contradiction. (In 
fact the inequality is strict, J2{(ij)\i^x*,jeY*} c ij > 0^ because f max > 0, and c^- and Cji are 
simultaneously zero or strictly positive.) 

b) (s E X*, t E X*). In this case the flow conservation implies that the net flow into 
Y* must be zero, i.e. T.{( ij )eA\iex*,jeY*}(fij ~ fji) = 0. Therefore E{fe)eA|iex*,jey*} r ij = 
12{(ij)eA\iex*,j£Y*} c ij — an d J2{(ij)eA\iex*,jeY*} r ji — J2{(ij)eA\iex*,jeY*} c ji > 0- As a result 
3 (ij) E A with = and Cji > 0, which is a contradiction. The case c) is similar to b). 
From all the above it follows that s E X* and t E Y* so (X*, Y*) is a cut. 

Now we have to prove that (X*,Y*) is also a minimum cut. We have 
E{(ij)eA\iexjeY} r ij = and therefore 

C(X,Y)= Cy= E (h ~ fji) = fmax (19) 

{(ij)eA\iex,jeY} {{ij)eA\iex,jeY} 
by the conservation of the flow. This implies, according to the Ford- Fulker son theorem, that 
(X, Y) is a minimum cut. 

Proposition 3: At the end of the algorithm A sg = A mc . 

Proof: A mc C A sg follows from the construction of the algorithm and Propositions 1 
and 2. In order to prove A sg C A mc let us assume, without any loss of generality, that the 
supergraph has a single cluster (multiple clusters are independent of each other). Let (IJ) 
be an arc of Q(Af,A) connecting subclusters I and J. ( We only consider the non-trivial 
case I ^ Y t and J ^ I s ; if I = Y t or J = X s the result follows immediately from the 
construction of the S cut and T cut.) Let us define: 

y*j = {all subclusters that can be reached from J along a directed path in Q(Af,A)} (20) 
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Xi—M — Yj (21) 

Now / ^ y,j, otherwise G(N,A) would contain a directed cycle. Therefore I G Xj, so 
(<-fc/,X0 is a directed partition in Q{M,A) and (/J) connects this directed partition. This 
partition determines a directed partition in (7™ aa: and therefore a minimum cut in G(N, A) 
according to Proposition 2. Then (IJ) C ^4 mc and because (IJ) is arbitrary A sg C *4 mc 
follows. 

Proposition 4- A mc is the set of arcs of G(N, A) that will be saturated by all the maximal 
flows (the 'weak' links of the network). 

Proof. Let (ij) G A mc and Gf ax the residual graph for a maximal flow. By the definition 
of A mc there exists a minimum cut (X, Y) such that i G X and j e F. This minimum cut is 
a directed cut in G™^ according to Proposition 2, so = 0, i.e. the arc (ij) is saturated. 
Conversely, let the arc (ij) with capacity c^- be saturated by all the possible maximal flows 
through the network, that is flows with value / = f max . Let us assume now that (ij) is not 
part of any minimum cut. Therefore, if we decrease the capacity by a small value e > 0, 
Cij — > c*j = Cij — e > 0, the maximum flow through the network will still have value f max , 
according to the Ford- Fulker son theorem. However, none of the maximal flows will now 'fit' 
through the network, in particular through the arc (ij). This is a contradiction, therefore 
the arc (ij) must be part of a minimum cut, that is (ij) G A mc . 
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FIGURES 

FIG. 1. Simplified representation of the network after the first step of the algorithm, showing 
the S cut, the T cut and Z (we denote X s by S and Y t by T - see text). Note that a directed arc 
connecting S and T may or may not exist. 

FIG. 2. Simplified representation of the network after the clusters making up Z have been 
identified and an example of a minimum cut. 

FIG. 3. The supergraph G{N, A) at the end of the algorithm and an example of a minimum cut. 

FIG. 4. Ground state structure of the bimodal RFIM for h/J = 3/2. The spins frozen 'up' are 
green, the spins frozen 'down' are white, while the other colors represent the spins making up Z\ 
neighbouring subclusters have different colors (see text). Note that Z is made up of independent 
clusters that contain one or more subclusters. A dot indicates a field 'up' and the absence of it 
indicates a field 'down'. 

FIG. 5. Ground state structure of the DAFF for h/J = 7/2 and p = 0.9. The color coding is 
the same as for the RFIM, only now black indicates an empty ('nonmagnetic') site. Note that in 
this case 'flipping' a (sub) cluster means coloring it with one or the other of the checkered patterns. 

FIG. 6. A RFIM cluster from Fig. || and its subcluster graph representation; R stand for red, 
Y for yellow and B for blue. The associated degeneracy is 7. 

FIG. 7. Interface configuration for the BDIM at p = 0.64. The spins frozen 'up' are green, the 
spins frozen 'down' are gray and the pieces that are not part of the percolating cluster (see text) 
are white. The other colors represent again the spins making up Z, with neighboring subclusters 
having different colors. Note that for this particular sample Z is made up of six independent 
clusters. 

FIG. 8. A different sample than in Fig. 6, that has three independent interfaces, each with its 
own excitations. Note that Z is in this case a single cluster in the algorithm denomination. 
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